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The use of APE smearing or other blocking techniques in fermion actions can provide many advantages. There 
are many variants of these fat link actions in lattice QCD currently, such as FLIC fermions. The FLIC fermion 
formalism makes use of the APE blocking technique in combination with a projection of the blocked links back 
into the special unitary group. This reunitarisation is often performed using an iterative maximisation of a gauge 
invariant measure. This technique is not differentiable with respect to the gauge field and thus prevents the use of 
standard Hybrid Monte Carlo simulation algorithms. The use of an alternative projection technique circumvents 
this difficulty and allows the simulation of dynamical fat link fermions with standard HMC and its variants. 



1. INTRODUCTION 

Fat Link Irrelevant Clover (FLIC) fermions 
have shown a number of promising advantages 
over standard actions, including improved con- 
vergence properties jjj and 0(a) improved scaling 
without the need for nonperturbative tuning [2]. 
Furthermore, a reduced exceptional configuration 
problem has allowed efficient access to the light 
quark mass regime in the quenched approxima- 
tion [3] . However, recent progress in the field has 
shown that the behaviour of quenched QCD can 
differ from the true theory both qualitatively and 
quantitatively in the chiral regime |4l5j . Thus, 
interest has now focused on dynamical QCD, be 
it (truly) unquenched, or partially quenched. 

In particular, as it is in the chiral regime where 
the difference from the valence approximation will 
be highlighted, we would like to go to light quark 
masses in dynamical QCD. This is an extremely 
computationally expensive endeavour. One might 
hope that the excellent behaviour at light quark 
mass displayed by FLIC fermions will carry over 
from the quenched theory to the unquenched one. 
This brings us to the issue of generating dynam- 
ical gauge field configurations with the fermionic 
determinant being that of the FLIC action. 

2. MaxReTr PROJECTION 

The standard technique for simulating dynam- 
ical fermions has for some time now been Hy- 



brid Monte Carlo (HMC) 0. The HMC algo- 
rithm alternates between a global Metropolis ac- 
cept/reject step and a gauge field update through 
the use of a hybrid Molecular Dynamics integra- 
tion, 

U m (x,t+At) = U„(x,T)exp(iATP^(x,T)), (1) 
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Here, the P^ix) are the momenta canonically con- 
jugate to the gauge field. The update of the mo- 
menta involves the derivative of the action with 
respect to the gauge field, a quantity which is also 
made use of in efficient local update algorithms. 

When attempting to use HMC with FLIC 
fermions, or some other fat link action for that 
matter, it is the calculation of Jjj which has 
proven problematic. FLIC fermions are clover- 
improved fermions where the irrelevant operators 
are constructed using APE smeared links |7I8| . 
These fat links are constucted by performing sev- 
eral sweeps of APE smearing across the lattice. 
Each sweep consists of an APE blocking step, 
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Sweep 



Unit Circle 



MaxReTr 



0.866138301214314 

1 0.960313394813806 

2 0.980735000838119 

3 0.988384926461589 

4 0.992103013943516 

5 0.994182852413813 

6 0.995457365275018 

7 0.996293668622924 

8 0.996878305318083 



0.866138301214314 

0.960348747275940 
0.980751346847750 
0.988393707639555 
0.992107844842705 
0.994185532052157 
0.995458835653863 
0.996294454083006 
0.996878710433084 



Table 1 

The mean link uq for a single configuration as a function of number of APE smearing sweeps, for the two 
different projection methods. The boldface indicates significant digits which match. The configuration is 
a dynamical gauge field with DBW2 glue and FLIC sea fermions, at f3 = 8.0, k = 0.128. 



followed by a projection back into ST/ (3) of the 



links, ujrHx) = V{V^'{x)). Typically, this pro- 
jection is performed by iteratively maximising the 
following gauge invariant measure, 



r(n) 



U^ l \x)[U {n - 1) ] 



max ReTr(t/'(a;)V; / (n)t (a;)). 



(4) 



Any fermion action which makes use of this re- 
unitarisation method cannot be simulated using 
HMC, as the use of the iterative algorithm pre- 
vents one from being able to construct 



(breaking the Z3 ambiguity by choosing the prin- 
cipal value of the cube root). Using the mean link 
as a measure of the smoothness of the smeared 
gauge field, Table ^ indicates that the two meth- 
ods presented here are nearly equivalent. 

However, the matrix inverse square root func- 
tion can be approximated by a rational polyno- 
mial (whose poles lie on the imaginary axis), 



W[X] « W k [X] = d X(X^X+c )J2 



1=1 



At A - 
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3. UNIT CIRCLE PROJECTION 

Given any matrix A, then X'X is hermitian 
and may be diagonalised. Then it is possible (for 
det A 7^ 0) to define a matrix [5| 
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X 



VAtA 



(5) 



whose spectrum lies on the complex unit circle 
and is hence unitary. Furthermore, W possesses 
the same gauge transformation properties as A. 
We can then constuct another matrix, 



w = 



1 



VdetW 



-.W 



(G) 



which is special unitary. This technique pro- 
vides an alternative method of SC/(3) projection 



This approximation is differentiablc in a matrix 
sense for all A for which the inverse square root 
can be defined. This means that we can construct 
J|j for fermion actions which involve unit circle 
projection. 

4. RESULTS 

Having now written down the APE smearing 
prescription (with projection) in a differentiablc 
closed form, the equations of motion necessary 
for the use of the HMC algorithm can be de- 
rived. Even in the simple case of APE smearing 
these equations are quite complex, but we will say 
something here about the mechanism for deriving 
them. Given the FLIC fermion action Sflic(^0> 
which has explicit dependence upon the thin links 
U and the fat links U', it is straightforward to 
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Sgaugc 


Pace 


a 




4.2 


0.1246 


IMP 


0.86 


0.112 


0.813 


4.3 


0.1253 


IMP 


0.82 


0.096 


0.727 


4.4 


0.1255 


IMP 


0.88 


0.092 


0.616 


4.5 


0.1253 


IMP 


0.82 


0.080 


0.631 


4.6 


0.1254 


IMP 


0.83 


0.077 


0.591 


9.0 


0.1224 


DBW2 


0.78 


0.144 


0.982 


9.5 


0.1228 


DBW2 


0.82 


0.114 


0.872 


10.0 


0.1234 


DBW2 


0.82 


0.105 


0.755 


10.5 


0.1236 


DBW2 


0.78 


0.099 


0.740 


11.0 


0.1239 


DBW2 


0.80 


0.092 


0.643 



Table 2 

The acceptance rate, lattice spacing and pion 
mass (in lattice units) for different dynamical sim- 
ulations, using (two flavours of) FLIC sea quarks 
and both Lusher- Wiesz (IMP) glue and DBW2 
glue. These preliminary results are obtained from 
20 12 3 x 24 configurations. Simulations are done 
using multiple time step HMC with trajectories 
of unit length, Arp SGudo f ormion = 0.02, Ar gaugc = 
0.01. 



write down the partial matrix derivatives ^ and 
■§lj7- However, to write down in one step the to- 
tal matrix derivative 4^ would be difficult, and 

all ' 

inefficent numerically Instead, one can use some 
calculus to both simplify the process and make 
it numerically efficient. This is done through the 
use of the matrix chain rule, 

dU ~ W + dU'*~dU' ( ' 

The full details of the derivation are presented 
elsewhere |10|. We find that the numerical over- 
head for the equations of motion are small, and 
the generation of dynamical gauge fields with 
FLIC fcrmions is still dominated by the neces- 
sary conjugate gradient inversion required as part 
of the pseudo-fermion formulation. Some prelim- 
inary results are presented in Table 

5. CONCLUSION 

Unit circle projection can be written in closed 
form and represented by a rational polynomial 
approximation. The matrix rational polynomial 



is differentiable with respect to the gauge field 
and thus enables one to derive the equations of 
motions necessary for the use of HMC with FLIC 
fermions. This provides a significant advantage 
over the MaxReTr method, where the absence of 
a differentiable form precludes the use of HMC 
and hence the availability of a simulation algo- 
rithm that scales linearly with the lattice volume 
V, although there are 0{V 2 ) alternatives [TT] , 
This means that we now have an 0(V) algorithm 
available for the dynamical simulation of FLIC 
fcrmions. Furthermore, the method is general 
and can be applied to any fermion action with 
reuniterisation, including overlap fermions with 
a fat link kernel |12I13I14I15| . or other types of 
fat link actions ^B] that may involve alternative 
smearing techniques |17| . 
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